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Abstract 



We report studies of the response of a massless particle confined by a potential. At large momen- 
. turn transfer q it exhibits y or equivalently Nachtmann ^ scaling, and acquires a constant width 

independent of q. This width has a contribution from the final state interactions of the struck 



Q ' particle, which persists in the |q| — > oo limit. The width of the response predicted using plane 

^ ■ wave impulse approximation is smaller because of the neglect of final state interactions in that 

- ^ , approximation. However, the exact response may be obtained by folding the approximate response 

X 

5h ' with a function representing final state interaction effects. We also study the response obtained 

from the momentum distribution assuming that the particle is on the energy shell both before 
and after being struck. Quantitative results are presented for the special case of a linear confining 
potential. In this case the response predicted with the on-shell approximation has correct values 
for the total strength, mean energy and width, however its shape is wrong. 
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I. INTRODUCTION 



Scattering of high energy probes from composite systems, such as electron scattering 
by nuclei [Q] and nucleons 0, or neutron scattering by liquid helium 0, is often used 
to study the structure of the bound system. The common assumption is that in deep 
inelastic scattering (DIS) at sufficiently high energy the probe is incoherently scattered by 
the constituents of the system. In the plane wave impulse approximation (PWIA), which 
neglects final state interaction (FSI) effects, DIS is directly related to the momentum and 
energy distribution of the constituents in the target. 

The role of FSI effects has been studied extensively in electron scattering from nuclear 
targets @, |^ and neutron scattering from liquid helium P]. Recently it has been suggested 
that they may also influence DIS of leptons by hadrons f^. In the present study we focus 
on scattering from targets with conflned constituents. The corresponding physical case 
concerns DIS from nucleons where, in distinction from the nuclear and liquid helium cases, 
the constituents are conflned in both the initial and flnal states. 

We have recently studied the response of a massless particle conflned by a linear 
(flux-tube) potential using the simple model Hamiltonian 

H=^/\^+V^r, (1) 

to gain insights into DIS of leptons by hadrons. The response to a hypothetical scalar probe 
is calculated as a function of the momentum and energy transfer, q, u in the lab frame. It 
is given by 

/?(q, u) = Y,\{I\e'''"'m'6{Ej -Eo- u), (2) 

where the sum is over all energy eigenstates /. The natural scaling variable in the many- 
body theory approach to DIS is y = z/ — |q| In the large |q| or scaling limit the response 
depends only on y, not on |q| and u independently. This scaling of the response is equivalent 
to the Nachtmann scaling, since the Nachtmann |0 scaling variable ^ = —y/M, where M 
is the hadron mass. In Ref. the -R(q, z/) is calculated for |q| < 10 GeV for the typical 
value of the string tension ^/o' = 1 GeV/fm in QCD, by calculating all the relevant states 
|/) contributing to the response. 

In this work we study the effects of the FSI of the struck particle on the response; analytic 
calculations of the width of the response are presented in Sec.|I|, and the numerical results 
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for a linear confining potential are given in Sec.|nT[ Both indicate that the FSI increase 
the width of the response beyond that predicted by PWIA. The analytic calculations also 
consider the nonrelativistic problem, in which q is large compared to all the momenta in 
the target, but smaller than the constituent mass m. The main differences between the 
nonrelativistic and the relativistic response are that the former peaks at z/ = |qp/2m and 
has a width proportional to |q|, while the latter peaks at z/ ~ |q|, and has a constant width 
in the scaling limit. The folding function representing the effects of the FSI on the 



response is discussed in the last section, Secl^, where conclusions are given. 



II. MOMENTS OF THE RESPONSE 

In the case of a single confined particle, the state of the system after the probe has struck 
the target is 

\X) = e^'i-^lO), (3) 

where |0) denotes the ground state of the particle. The state |X) is not an eigenstate of 
the Hamiltonian and therefore has a distribution in energy. It has a unit norm, {X\X) = 
(0|e~*'' '"e*^'''|0) = 1. The total strength of the response, given by the static structure function 

POO 

5(|q|)= / duRilqlu), (4) 
Jo 

is therefore unity. In many-body systems 5'(|q|) is not necessarily equal to one. Subsequent 
formulas pertain to the general case and show factors of 5'(|q|) explicitly. 

The mean excitation energy of the state |X) is given by the first moment of the response: 

= sMiXlH - E„\X) = ^ fd..R(H,.). (5) 

The width of the distribution in energy is characterized by the second moment of the energy 
about the mean: 



A^(|q|) = ^(^1 (h - ^^TT^^ ) 1^)- (6) 
For example, if the response is Gaussian, it is completely determined by these three 
moments: 

The full-width at half-maximum (FWHM) of a Gaussian response is given by2A(|q|)yM2 ~ 
|A(|q|). Obviously many more moments are required to describe a response of more general 
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shape. Nevertheless, a necessary condition for the occurrence of the y or Nachtmann ^ scahng 
is that 5'(|q|), z^(|q|) — |q|, A(|q|), and all the higher moments become independent of |q|, 
since the response depends only on v — |q| when |q| oo. 

In the remainder of this section we will calculate the average excitation energy, X^dq]) 
and the width, A(|q|) for the model Hamiltonian of Eq.(|l|) in the scaling limit, |q| — > oo. 
They are calculated analytically for a general potential V{r). Section |T| presents numerical 
results for the case of a linear confining potential. This analysis allows the comparison of 
the results for the exact response with the following two common approximations. The 
PWIA, which assumes that the final state of the struck particle can be approximated by a 
plane wave in an average potential, is commonly used in nuclear physics. And the on-shell 
approximation (OSA), in which the struck particle is assumed to be in a plane wave state 
with the energy of that of a free particle, before and after the interaction with the probe, 
used in hadron physics. 

We proceed with the evaluation of the average excitation energy [Eq.(|)] for the case of 
the exact response. The kinetic energy term of the matrix element {X\H\X) is 

^rz(A:)|k + q|, (8) 

where we have transformed to momentum space and n{k) = \ipo{k)\'^. In the scaling limit, 
with q = |q|z, we can expand |q + k| and obtain: 

|k + q| = |q| + k, + ^^kl - ^^klk, + O . (9) 

Here O (^j^j denotes the neglected terms of that and higher order, and k\ = + ky. 
Substituting this expansion into Eq.@ yields, 

^(|q|) = |q| + (^)o -Eo + ^(^')o + O ; (10) 

the term of O (^j^^ is zero by symmetry. Note that Eq — {V)o = (T)o, where T denotes 
the kinetic energy. Thus I7(|q|) = |q| — {T)o in the limit |q| oo. The requirement that 
'^(IqI) ~ |q| becomes constant is naturally satisfied in this limit. 

The width [Eq.(^] depends on the matrix element {X\H'^\X)] = + \/^V{r) + 
V{r)^/j^ + V'^{r). The first term of gives: 

(X|/|X) = (|k + q|2)o = |qp + (A;2)o. (11) 
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The average of the second and third terms in H"^ may be shown to be equal to each other 
and 

^ '^'^^^°+2^^^'^^° + ^(l^)' ^^^^ 
where [K-?/'o](^) = / d?'re^^'^V{r)%ljQ{r). Substitution into Eq.(|^) with 5'(|q|) = 1 gives: 

A^(|q|) = + {V\ - {V)l + ^ ((fcV)o - {k\{V),) + O . (13) 

This expression demonstrates that the width of the exact response is independent of |q| 
in the hmit |q| — ^> oo as necessary for y scahng. It also shows that the width has a kine- 
matic contribution dependent upon the target momentum distribution, and an additional 
interaction contribution. 

As mentioned, the PWIA assumes that a constituent of momentum k, after being struck 
by the probe, may be described by a plane wave with momentum k + q in an assumed 
average potential chosen to give the exact V of Eq.(p!oD. The PWIA response is 

^n(A;)5(e(k + c^) - E, - u) (14) 

where e(k + q) is the energy of the plane wave, taken to be: 

e(k + q) = |k + q| + (y)o. (15) 

The first term is the kinetic energy of the struck particle. The second term is the average of 
the potential in the final state of the system given by, (X|1^|X) = (0|e^*'^'''V^e*'' '"|0) = {V)q. 
The Rpwia{\^\iV) exhibits y scaling as can be easily seen by expanding the argument of 
the 5 function in RpwiA- In the large |q| limit RpwiA depends only on y = — |q|. 
The average energy and the width of the PWIA response are calculated using: 

dvv^RpwiA{H,i^) = duu^ J -^n{k)6{eik + q)-Eo~u). (16) 

The average excitation energy, obtained with n = 1, agrees with the exact result in Eq.(p!OD 
by construction. 

The width of the PWIA response, however, is not the same as the width of the exact 
response. 

Al^j^ = 1(^)0 + O(r^) (17) 
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contains only the first term of tlie exact result [Eq.(0)] due to the target momentum dis- 
tribution. The second term, {V^)o — {V)l, of the exact represents the FSI contribution 
neglected in the PWIA. It does not vanish in the |q| ^ oo limit for relativistic kinematics. 

I |2 

In the non-relativistic case, H^ji = ^ + V{r), the exact V is given by: 

= i^lqP . (18) 
The nonrelativistic-PWIA (NR-PWIA) response: 

RNR-PwiA{\q\,iy) = J -0y^^{k)6 (^-^^^ + {V)o -Eo-u^ (19) 
gives the exact V. For the width of the NR-PWIA response we obtain: 



Note that in Eqs.(|T^) and (pO|) we have not taken the |q| — >■ oo limit. 
The width of the exact NR response is: 

A^^(|q|) = ^In^PwiA + ^ {{vk\ - iyuk\) + iy\ - iy)l . (21) 

It differs from ^nr^pwia in terms of order l/|q| which can be neglected in the scaling limit. 
Thus, in contrast to the relativistic case, the FSI do not increase the width of the NR-PWIA 
response at large |q|. 

Finally we consider the on-shell approximation (OSA) in which the energy of the struck 
constituent is that of a free relativistic particle before and after the interaction with probe, 
as assumed in the quark-parton model. The response in OSA is 

^n(A;)5(|k + q|-|k|-z/); (22) 

it depends only on the momentum distribution of target constituents and obeys y scaling. 
The average excitation in OSA is 

VosA = |q| - (T)o + :^{k')o + O , (23) 

3|q| Vlql/ 

and the width is given by: 

Ao5A = + {k\ - {T)l + o(^-^y (24) 
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FIG. 1: The response versus y calculated exactly for Tq = 100 MeV (thin solid curve) and Tq = 50 
MeV (dotted curve). The response in OS A are shown for |q| = 10 GeV (dashed) and |q| ^ oo 
(dot-dashed). The PWIA response for |q| = 10 GeV and |q| ^ (X) lie on essentially the same (thick 
solid) curve. 

The exact value of V [Eq.(|TO|)] is reproduced by the OSA for any potential. However, the 
^osA (fc^)o — (T)q in place of the {V^)o — {V)l in the leading term of the exact 
[Eq.(|T3|)]. For a massless particle in a linear confining potential, i.e. for the Hamiltonian of 
Eq.(|T]), (T)o = {V)o, and {k'^)o = (^^)o- Therefore for this particular Hamiltonian the OSA 
reproduces the exact value of A. 

The calculated responses shown in the next section however indicate that the shape of 
the OSA response is incorrect. We have calculated the third moment of the H in the state 
|X) and observed that the exact and OSA responses indeed give different values. We expect 
the moments to differ also for orders higher than the third since the shapes of the exact and 
the OSA response are very different. 

III. NUMERICAL RESULTS 

We first compare the response functions for |q| = 10 GeV before comparing their mo- 
ments. In Ref. it has been shown that the scaling limit is obtained for such values of 
|q|. The exact response, Eq.(|[), is a sequence of 6 functions at u = Ej — Eq. In order to 
obtain a smooth response we assume decay widths Fq for all the excited states. Note that 
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the energies of the states |/) that contribute to the response at |q| = 10 GeV are large, 
therefore their decay widths are not affected by the energy dependent terms assumed in Ref. 
1^. The response including decay widths is given by: 

The responses obtained with Tq = 100 and 50 MeV are shown in Fig. ^ along with the 
PWIA and OS A responses for |q| = 10 GeV and for |q| oo. The difference between the 
exact responses for Fq = 100 and 50 MeV are much smaller than those between the exact 
and the approximate. 

We note that the shape of the PWIA response is qualitatively similar to that of the exact, 
however, its width is too small. This is a direct consequence of the neglect of interaction 
terms in A [Eq.(|l3|)] as discussed in the last section. The width A of the response for Fq 0, 
is 409 MeV, while the Ap^^iA = 326 MeV. Note that for a Gaussian response, Eq.(^ the 
FWHM is ~ 5 A/3. The FWHM of the exact and the PWIA responses shown in Fig. |] are 
larger because they are not exactly Gaussian. 

The OSA results in the discontinuous curves shown in Fig. |l[ They are discontinuous 
at the lightline (|q| = u) because the response of free particles is limited to the spacelike 
region u < |q|. The discontinuity at y = is in clear conflict with the exact response which 
is continuous across the lightline and is non-zero in the timelike {y > 0) region. Therefore 
the OSA appears to be unsatisfactory even though for the special case of a linear potential 
it has the exact values of 5'(|q|), I^(|q|) and A(|q|). 



IV. FSI FOLDING FUNCTION 



The exact response may be obtained from Rpwjj^[\q\, u) via convolution with the folding 
function 0, ^, The states |k + q) used in the PWIA are not eigenstates of the H. 
Therefore they must have a distribution in energy, and we may express the exact response 
as: 

R{\q\,u) = j dv'F{\^lv - v')Rpwia{H,i^'). (26) 

The folding function F defined above, is primarily meant to describe the width of the plane 
wave states. Scaling occurs even in presence of interactions because this folding function 



becomes independent of |q| at large |q| |10, 11]. Occurrence of scaling does not imply that 
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FIG. 2: Folding functions for large |q|. The solid and dashed curves correspond to Fq = 50 and 



the folding functions obtained from the exact and PWIA response by deconvolution are peaked at 
V ~ 0.2 GeV. 

either the PWIA or the OSA is valid. Fig. |^ shows the folding function obtained from the 
PWIA and exact responses for |q| = 10 GeV and Fq = 100 and 50 MeV. They are calculated 
numerically from -R(|q|, v) and -RpvK/yi(|q|) ^) by deconvolution. 

In order to extract information on the ground state wave function from the measured 
response, we must first remove the FSI effects represented by the folding function F(|q|, z/). 
However, in most cases the folding function is not known. In some nonrelativistic cases it 
can be calculated with quantum Monte Carlo methods [T^. It can also be estimated using 
the Glauber approximation [|], [T^. Here we consider the possibility of approximating the 
F by a Gaussian folded with the Lorentzian distribution used in Eg . (|25|) . In this case the 
width of the Gaussian is given by the interaction contribution to the exact A Eg . (p!3D : 



100 MeV, respectively. The approximate folding functions Fa of Eq.(p9|) are peaked at = while 




(27) 
(28) 



and the convolution of Fq with the Lorentzian is given by: 




(29) 



(30) 
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FIG. 3: The exact response (solid curve) compared with the response obtained by folding the 
PWIA response with the approximate folding function Fa of Eq.(|2^) for Tq = 100 MeV. 

The Fa is shown in Fig. 0. It does not have the finer structures in the exact F, and it peaks 
at z/ — z/' = 0, while the F peaks at ~ 0.2 GeV. Nevertheless, a much better approximation 
to the response is obtained by folding the PWIA response with the Fa as shown in Fig. ^ 
In conclusion, we find that the effects of FSI, neglected in the PWIA, become independent 
of |q| at large |q| for a particle confined in a linear potential. They increase the width of the 
response, and must be removed before extracting momentum distributions or other target 
structure information from the observed response. In general one can always represent the 
FSI effects by a folding function; in this simple case a Gaussian folding function provides 
a reasonable approximation. Finally, the OSA predicts a response of the wrong shape but 
with correct values for the first three moments in the case of a linear potential. 
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